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ABSTRACT 


Independent Earth orientation measurements taken by the space-geodetic techniques of 
lunar and satellite laser ranging, very long baseline interferometry, and the Global Positioning 
System have been combined using a Kalman filter. The resulting combined Earth orientation 
series, SPACE2013, consists of values and uncertainties for Universal Time, polar motion, and 
their rates that span from September 28, 1976, to June 30, 2014, at daily intervals and is available 
in versions with epochs given at either midnight or noon. The space-geodetic measurements used 
to generate SPACE2013 have then been combined with optical astrometric measurements to form 
two additional combined Earth orientation series: (1) COMB2013, consisting of values and 
uncertainties for Universal Time, polar motion, and their rates that span from January 20, 1962, 
to June 30, 2014, at daily intervals and which are also available in versions with epochs given at 
either midnight or noon; and (2) POLE2013, consisting of values and uncertainties for polar 
motion and its rate that span from January 20, 1900, to June 22, 2014, at 30.4375-day intervals. 
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INTRODUCTION 


Reference series of Earth orientation parameters (EOPs) obtained by combining 
independent measurements of the Earth’s orientation are generated annually at the Jet Propulsion 
Laboratory (JPL) in support of tracking and navigation of interplanetary spacecraft. This report 
describes the generation of the most recent such combined Earth orientation series: SPACE2013, 
COMB2013, and POLE2013. Since the procedures used to generate these most recent series are 
similar to those used to generate previous such combinations, only a brief description of their 
generation is given here. Further details regarding the approach used at JPL to annually combine 
independent measurements of the Earth’s orientation can be found in Gross ( 1 996, 2000) and Gross 
et al. (1998). 
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SPACE2013 


Data Sets Combined to Form SPACE2013 

SPACE2013 is a combination of independent space-geodetic measurements of the Earth’s 
orientation. Table 1 lists the space-geodetic series used in generating SPACE2013, giving their 
identifiers, the number of measurements from each series that were actually incorporated into 
SPACE2013, and the time interval spanned by those measurements. Note that the University of 
Texas Center for Space Research (UTCSR) satellite laser ranging (SLR) Universal Time (UT) 
values were not used in generating SPACE2013 due to problems associated with separating this 
component of the Earth’s orientation from the effects of unmodeled forces acting on the satellite 
that cause the node of its orbit to drift (see Gross et al. 1998, p. 217 for further discussion about 
this point). For similar reasons, the International Laser Ranging Service (ILRS) satellite laser 
ranging length-of-day (LOD) values have also not been used in generating SPACE2013. 

Since it was desirable to combine only independent measurements of the Earth’s 
orientation, only those series listed in Table 1 were used, even though other space-geodetic series 
are available from other analysis centers. When more than one series determined by the same 
measurement technique was used, care was taken to make sure that the measurements themselves 
were not included more than once. In particular, polar motion measurements from the JPL Global 
Positioning System (GPS) series were only used until the start of the International Global 
Navigation Satellite Systems (GNSS) Service (IGS) combined series EOP(IGS) 95 P 01 on 
January 1, 1995; polar motion measurements from the IGS combined series EOP(IGS) 95 P 01 
were then used until the start of the IGS combined series EOP(IGS) 95 P 02 on June 30, 1996; 
polar motion measurements from the IGS combined series EOP(IGS) 95 P 02 were then used until 
the start of the accumulated IGS Solution Independent Exchange (SINEX) combined series 
EOP(IGS) 00 P 03 on February 26, 2000; and polar motion measurements from the accumulated 
IGS SINEX combined series EOP(IGS) 00 P 03 were used thereafter. Similarly, LOD 
measurements from the IGS combined series EOP(IGS) 95 P 02 were used until it ended on 
February 25, 2000 with LOD measurements from the accumulated IGS SINEX combined series 
EOP(IGS) 00 P 03 used thereafter. In addition, UT measurements from the National Oceanic and 
Atmospheric Administration (NOAA) International Radio Interferometric Surveying (IRIS) 
Intensive UT1 series were used until it ended on December 31, 1994; measurements from the 
United States Naval Observatory (USNO) National Earth Orientation Service (NEOS) Intensive 
UT1 series were then used until it ended on December 4, 2000; and measurements from the 
Goddard Spaceflight Center (GSFC) NEOS Intensive UT1 series were used thereafter. Finally, 
polar motion measurements from the UTCSR SLR series EOP(CSR) 96 L 01 were used until it 
ended on February 3, 1996 with measurements from the ILRS combined SLR series being used 
thereafter. 

Data Preprocessing and Treatment of Tide-Induced Rotational Variations 

The Earth orientation series listed in Table 1 were first preprocessed by removing leap 
seconds from the Universal Time (UT1) values and, when necessary, by correcting the UT1 values 
to be consistent with the extended definition of Greenwich Sidereal Time (GST) as adopted by the 
International Earth Rotation and Reference Systems Service (IERS; IERS 1997, p. 149). Since 
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most of the series listed in Table 1 were already consistent with the extended definition of GST, 
this correction needed to be applied to only the NOAA IRIS Intensive UT1 series. 


Table 1 . Data Sets Combined to Form SPACE2013* 


Data Set Name 

Data Type 

Analysis 

Center 

Reference 

Data Span 

Number 

LLR (JPL14M01; VOL, UTO) 
McDonald Cluster 

LLR 

JPL 

Williams et al. (2012) 

Oct. 5, 1976, to Jan. 5,2012 

747 

OCA 

LLR 

JPL 

Williams et al. (2012) 

April 7, 1984, to June 8, 2014 

1218 

Haleakala 

LLR 

JPL 

Williams et al. (2012) 

Feb. 10, 1985, to Aug. 11, 1990 

70 

Apache Point 

LLR 

JPL 

Williams et al. (2012) 

June 4, 2006, to Sep. 1, 2013 

416 

UTCSR (CSR96L01) 
LAGEOS Polar Motion 

SLR 

UTCSR 

Eanes & Watkins ( 1996) 

Sep. 27, 1976, to Feb. 3, 1996 

2219 

ILRSA (29JUL14; Polar motion) 
ILRS Primary Combination 

SLR 

ASI-CGS 

Sciarretta et al. (2010) 

Feb. 4, 1996, to June 29, 2014 

6683 

DSN (JPL97R0 1 ; T, V) 
Califomia-Spain Cluster 

VLBI 

JPL 

Steppe et al. (1997) 

Nov. 26, 1979, to June 29, 2014 

1080 

Califomia-Australia Cluster 

VLBI 

JPL 

Steppe et al. (1997) 

Oct. 28, 1978, to Sep. 30, 1997 

698 

NOAA (NOAA95R02) 
IRIS Intensive UT1 

VLBI 

NOAA 

Ray et al. (1995) 

April 2, 1984, to Dec. 31, 1994 

2395 

USNO (N9903) 

NEOS Intensive UT1 

VLBI 

USNO 

Eubanks et al. ( 1999) 

Jan. 4, 1995, to Dec. 4. 2000 

1497 

NASA/GSFC (GSF2014a) 
NEOS Intensive UT1 

VLBI 

GSFC 

NASA/GSFC VLBI Group (2014) 

Dec. 6, 2000, to June 30, 2014 

4230 

NASA/GSFC (GSF2014a) 
Multibaseline 

VLBI 

GSFC 

NASA/GSFC VLBI Group (2014) 

April 12, 1980, to June 26, 2014 

4737 

NASA/GSFC (GSFC 1 122) 
Westford-Fort Davis 

VLBI 

GSFC 

Gordon et al. (1999) 

June 25, 1981, to Jan. 1, 1984 

105 

Westford-Moj ave 

VLBI 

GSFC 

Gordon et al. (1999) 

March 21, 1985, to Aug. 6, 1990 

18 

GPS (21APR04; Polar motion) 
Post-processed Flinn Analysis 

GPS 

JPL 

Fleflin et al. (2004) 

June 10, 1992, to Dec. 31, 1994 

817 

GPS (IGS95P01; Polar motion) 
IGS Final Combined (ITRF93) 

GPS 

NRCan 

Kouba & Mireault (1997) 

Jan. 1, 1995, to June 29, 1996 

546 

GPS (IGS95P02; Polar motion) 
IGS Final Combined (IGS08) 

GPS 

CODE 

Mireault & Kouba (2000) 

June 30, 1996, to Feb. 25, 2000 

1337 

GPS (IGS00P03; Polar motion) 
IGS SINEX Combined (IGS08) 

GPS 

IGN 

Rebischung & Garayt (2013) 

Feb. 26, 2000, to June 30, 2014 

5227 

GPS (IGS95P02; LOD) 

IGS Final Combined (IGS08) 

GPS 

CODE 

Mireault & Kouba (2000) 

Feb. 23, 1997, to Feb. 25, 2000 

1098 

GPS (IGS00P03; LOD) 

IGS SINEX Combined (IGS08) 

GPS 

IGN 

Rebischung & Garayt (2013) 

Feb. 26, 2000, to June 30, 2014 

5226 


* LLR, lunar laser ranging; JPL, Jet Propulsion Laboratory; VOL, variation of latitude; UT, Universal Time; OCA, Observatoire de la Cote d’Azur; UTCSR, 
University of Texas Center for Space Research; LAGEOS, Laser Geodynamic Satellites; SLR, satellite laser ranging; ILRS, International Laser Ranging 
Service; ASI, Agenzia Spaziale Italiana; CGS, Centro di Geodesia Spaziale; DSN, Deep Space Network; T, transverse; V, vertical; VLBI, very long baseline 
interferometry; NOAA, National Oceanic and Atmospheric Administration; ERIS, International Radio Interferometric Surveying; USNO, United States 
Naval Observatory; NEOS, National Earth Orientation Service; NASA, National Aeronautics and Space Administration; GSFC, Goddard Space Flight 
Center; GPS, Global Positioning System; IGS, International Global Navigation Satellite System (GNSS) Service; ITRF, International Terrestrial Reference 
Frame; NRCan, Natural Resources Canada; CODE, Center for Orbit Determination in Europe; IGN, Institut National de l'lnformation Geographique et 
Forestiere; SINEX, solution independent exchange; LOD, length of day. 
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Rotational variations caused by solid Earth and ocean tides were also removed from the 
UT1 values. The effect of the solid-Earth tides was removed by using the model of Yoder et al. 
(1981); the model of Kantha et al. (1998) was used to remove the effect upon UT1 of the ocean 
tides at the Mf, Mf, and Mm tidal frequencies. Since the Yoder et al. (1981) model already 
includes a contribution from the equilibrium ocean tides, only the Kantha et al. (1998) oceanic 
corrections to the Yoder et al. (1981) model were actually removed. Also note that the Kantha et 
al. (1998) model was used to remove the effect of ocean tides on only UT1, not on polar motion. 
Ocean-tide-induced polar motion variations were not removed from any of the polar motion 
observations. Finally, the only Earth orientation series listed in Table 1 that includes the effects 
of semidiurnal and diurnal ocean tides on the Earth’s orientation is the NOAA IRIS Intensive UT1 
series. This series included these effects by adding to the released UT 1 values the model of Herring 
(1993; also see Herring and Dong 1994). Hence, the same Herring (1993) model was used to 
remove them. 

On June 30, 1996, the IGS reference frame changed from the International Terrestrial 
Reference Frame ITRF93 to ITRF94; on March 1, 1998, it changed from ITRF94 to ITRF96; on 
August 1, 1999, it changed from ITRF96 to ITRF97; on February 27, 2000, it changed from 
ITRF97 to the IGS realization of ITRF97 known as IGS97; on December 2, 2001, it changed from 
IGS97 to the IGS realization of ITRF2000 known as IGS00; on January 11, 2004, it changed from 
IGS00 to a new IGS realization of ITRF2000 known as IGbOO; on November 5, 2006 it changed 
from IGbOO to the IGS realization of ITRF2005 known as IGS05 (Ferland, 2006); on April 17, 
2011 it changed from IGS05 to the IGS realization of ITRF2008 known as IGS08 (Rebischung et 
al., 2011, 2012); and on October 7, 2012 IGS08 was replaced by the new IGS realization of 
ITRF2008 known as IGb08 (Rebischung, 2012). These changes in reference frames potentially 
introduce discontinuities into the IGS combined Earth orientation series. However, the IGS 
Final combined series EOP(IGS) 95 P 02 and the accumulated IGS SINEX combined series 
EOP(IGS) 00 P 03 used here have had each of these segments aligned to the same IGS reference 
frame. Thus, to within the uncertainty in determining the corrections required to align each 
segment, the IGS combined series EOP(IGS) 95 P 02 and EOP(IGS) 00 P 03 used here should be 
reasonably consistent with each other. They were, therefore, concatenated with one common set 
of bias-rate corrections being determined for them, as described below. 

Adjustments Made to Space-Geodetic Series Prior to Combination 

Prior to combining the series listed in Table 1 to form SPACE2013, series-specific 
corrections were applied for bias and rate, the stated uncertainties of the measurements were 
adjusted by multiplying them by series-specific scale factors, and outlying data points were 
deleted. Values for the bias-rate corrections and uncertainty scale factors were detennined by an 
iterative procedure wherein each series was compared, in turn, to a combination of all others. In 
order to minimize interpolation error (see Gross et al. 1998, pp. 223-225), the comparison of each 
series to its reference combination was done at the epochs of the measurements of that series by 
generating its reference combination using a Kalman filter that interpolates to (and prints its EOP 
estimates at) the exact epochs of those measurements. Also, both the bias-rate corrections and the 
uncertainty scale factors for all components of a given series were determined simultaneously in a 
multivariate approach using nonlinear weighted least squares. Using a multivariate approach 
allows the correlations between the components to be taken into account when detennining the 
bias-rate corrections and uncertainty scale factors (see Gross et al. 1998, pp. 225). 
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All the series listed in Table 1 were included in the iterative procedure. Details of the 
iterative procedure are described in Gross (1996, 2000) and Gross et al. (1998) and will not be 
repeated here. These details include: 

(1) Use of a reference series, SPACE2012 (Ratcliff and Gross, 2013), for initial bias-rate 
alignment; 

(2) Analysis of each data type in its natural reference frame; 

(3) Clustering of the McDonald LLR stations and, separately, the DSN very long baseline 
interferometry (VLBI) stations in California, Australia, and Spain; 

(4) Initial convergence on values for the series-specific bias-rate corrections and uncertainty 
scale factors prior to the start of outlier detection and deletion; and 

(5) Final convergence on these values after detecting and deleting all data outliers. 

At the end of the iterative procedure, relative bias-rate corrections have been detennined 
that make the series agree with each other in bias and rate; uncertainty scale factors have been 
determined that make the residual of each series have a reduced chi-square near one when 
differenced with a combination of all others; and outlying data points (those with residual values 
greater than four times their adjusted uncertainties) have been deleted. A total of 105 data points, 
or about 0.26% of those combined, were thus deleted. 

For GSFC NEOS Intensive and GSFC Multibaseline VLBI, the initial alignment to 
SPACE2012 was performed individually for each baseline containing at least 90 data points and 
spanning more than 1 year. Baselines containing less than 1 year of data or fewer than 90 data 
points were aligned as a group to SPACE2012. Subsequent steps of the iterative procedure were 
then conducted independently on the individual baseline groups. 

Starting April 17, 201 1 the IGS changed their procedure for estimating uncertainty in the 
GPS EOP values. This introduced a change in the magnitude of the polar motion X (PMX), polar 
motion Y (PMY), and LOD uncertainties of approximately a factor of 6. To account for this sudden 
change, uncertainty scale factor values after April 17, 201 1 were estimated separately from those 
of earlier epochs using a trial reference series. These trial scale factors were applied to the original 
IGS GPS data which were then used in the iterative procedure described above to form a final 
reference series. The resulting scale factors shown in Table 2 for the IGS00P03 data set include 
the contribution of both the trial scale factors and the overall iteratively determined scale factor. 

Finally, each series was placed within a particular IERS reference frame by first 
applying to it an additional bias-rate correction to each series individually and then one final bias- 
rate correction that is common to all the series. Each series was compared to the IERS combined 
Earth orientation series EOP(IERS) 08 C 04 (Bizouard and Gambis, 2011) during the interval 
1984-2013 in order to obtain the additional corrections required to make it agree in bias and rate 
with the IERS series. The relative bias-rate corrections and uncertainty scale factors (determined 
by the iterative procedure described above) and the additional bias-rate corrections that align each 
series with the IERS series were applied to each individual series. These series were then 
combined and the resulting combination was also compared to EOP(IERS) 08 C 04 to obtain final 
bias-rate corrections that make the combined series agree in bias and rate with the IERS series. 
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Table 2. Adjustments to Space-Geodetic Series 1 


Data Set Name 


Bias. 

, mas 

Rate, mas/yr 

Uncertainty Scale Factor 

LLR (JPL14M01) 


VOL 

UT0 

VOL UT0 

VOL 


UT0 

McDonald Cluster 

-0.029 ± 0.086 

0.180 ±0.088 

0.029 ±0.020 -0.034 ±0.018 

1.478 


1.362 

OCA 


0.008 ± 0.035 

-0.002 ± 0.030 

-0.003 ± 0.004 0.00 1 ± 0.003 

2.240 


1.856 

Haleakala 


-2.052 ± 1.203 

-1.415 ±0.739 

-0.388 ±0.234 -0.271 ±0.155 

2.787 


2.421 

Apache Point 


-0.443 ± 0.384 

-0.668 ± 0.403 

0.020 ± 0.023 0.035 ± 0.024 

2.095 


2.142 

DSN (JPL97R01) 


T 

V 

T 

V 

T 


V 

Califomia-Spain Cluster 

-0.810 ±0.024 

-0.325 ± 0.059 

0.186 ±0.002 0.176 ±0.004 

2.099 


1.489 

Califomia-Australia Cluster 

0.607 ±0.018 

0.486 ±0.053 

0.078 ± 0.008 0.000 ± 0.024 

1.436 


1.168 

NASA GSFC (1122) 

T 

V 

T 

V 

T 


V 

Westford-Fort Davis 

8.559 ±3.838 

2.794 ±5.920 

0.826 ±0.386 0.240 ± 0.593 

1.247 


0.911 

W estford-Moj ave 


0.633 ± 1.097 

0.178 ± 1.673 

0.000 0.000 

2.500 


0.719 

GSFC (GSF2014a) b 

PMX 

PMY 

UT1 

PMX PMY 

UT1 

PMX 

PMY 

UT1 

G3GcKuRi 

0.006 ±0.102 

0.54 ±0.104 

0.227 ±0.168 

-0.06 1 ± 0.062 0.050 ± 0.064 

0.007 ±0.102 

1.346 

1.377 

1.899 

HrRiWfWz 

0.616 ±0.172 

-0.495 ±0.178 

0.127 ±0.284 

0.050 ±0.029 -0.1 16 ±0.030 

-0.002 ± 0.048 

1.174 

1.247 

1.701 

MoRiWfWz 

0.136 ±0.121 

-0.450 ±0.1 15 

-0.1 12 ±0.144 

-0.080 ± 0.052 -0.203 ± 0.048 

-0.037 ± 0.066 

1.448 

1.344 

0.376 

Other networks 

0.019 ±0.006 

-0.037 ± 0.005 

0.127 ± 0.010 

- 0.001 ± 0.000 - 0.001 ± 0.000 

-0.007 ±0.001 

1.577 

1.574 

1.642 

NOAA (95R02) 

PMX 

PMY 

UT1 

PMX PMY 

UT1 

PMX 

PMY 

UT1 

IRIS Intensive 



0.195 ±0.023 


0.038 ±0.007 



1.028 

USNO (N9903) 

PMX 

PMY 

UT1 

PMX PMY 

UT1 

PMX 

PMY 

UT1 

NEOS Intensive 



-0.129 ±0.040 


0.046 ± 0.007 



1.358 

GSFC (GSF2014a) b 

PMX 

PMY 

UT1 

PMX PMY 

UT1 

PMX 

PMY 

UT1 

KkNy 



-0.924 ± 0.554 


0.050 ±0.029 



1.206 

KkWz 



0.095 ± 0.022 


-0.005 ±0.001 



1.519 

TsWz 



0.202 ± 0.062 


-0.009 ± 0.004 



1.711 

NyTsWz 



0.019 ±0.215 


0.001 ±0.012 



1.832 

KkSvWz 



0.324 ±0.149 


-0.018 ±0.009 



1.754 

Other networks 



-0.310 ±0.022 


0.015 ±0.002 



1.598 

UTCSR (96L01) 

PMX 

PMY 

UT1 

PMX PMY 

UT1 

PMX 

PMY 

UT1 

LAGEOS 

-0.405 ±0.010 

0.734 ±0.009 


0.107 ±0.004 0.136 ±0.003 


0.765 

0.714 


ILRSA (29JUL14) 

PMX 

PMY 

LOD 

PMX PMY 

LOD 

PMX 

PMY 

LOD 

Primary Comb. 

0.053 ± 0.006 

-0.004 ±0.005 


-0.003 ± 0.000 -0.006 ± 0.000 


3.109 

3.005 


GPS (21APR04) 

PMX 

PMY 

LOD 

PMX PMY 

LOD 

PMX 

PMY 

LOD 

JPL Post-Flinn - 

-0.278 ±0.016 

0.141 ±0.013 


0.021 ±0.015 -0.032 ±0.012 


1.300 

0.954 


GPS (IGS95P01) 

PMX 

PMY 

LOD 

PMX PMY 

LOD 

PMX 

PMY 

LOD 

Final Combined - 

-0.146 ±0.056 

0.204 ± 0.054 


0.219 ±0.019 0.314 ±0.019 


1.300 

1.029 


GPS (IGS95P02) 

PMX 

PMY 

LOD 

PMX PMY 

LOD 

PMX 

PMY 

LOD 

Final Combined 

0.095 ± 0.004 

0.009 ± 0.004 

-0.019 ±0.010 

-0.005 ± 0.000 -0.001 ± 0.000 

0.002 ±0.001 

2.042 

1.635 

1.818 

GPS (IGS00P03) 

PMX 

PMY 

LOD 

PMX PMY 

LOD 

PMX 

PMY 

LOD 

before 4/17/2011 

0.095 ±0.004 

0.009 ±0.004 

-0.019 ±0.010 

-0.005 ± 0.000 -0.001 ± 0.000 

0.002 ±0.001 

2.042 

1.635 

1.818 

after 4/17/2011 

0.095 ± 0.004 

0.009 ± 0.004 

-0.019 ±0.010 

-0.005 ± 0.000 -0.001 ± 0.000 

0.002 ±0.001 

11.235 

5.960 

10.644 


a ) Reference date for bias-rate adjustment is 1993.0. See Table 1 footnotes, mas, milliarcseconds; PMX, polar motion X; PMY, polar motion Y. 

b) GSFC VLBI Station IDs from ftp://cddis.gsfc.nasa.gov/vlbi/ivscontrol/ns-codes.txt 

G3 = NRAO 85-3; Gc = Gilmore Creek; Hr = HRAS at Fort Davis; Kk = Kokee Park; Ku = 9m at Kokee Park; Mo = 12m at Goldstone; Ny = Ny 
Alesund, Svalbard, Norway; Ri = Richmond, FL; Sv = Svetloe, Russia; Ts = Tsukuba, Japan; Wf = Westford, MA; Wz = Wettzell, Germany. 
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The total bias-rate corrections and uncertainty scale factors determined by the procedures 
outlined above are given in Table 2. The values of the bias-rate corrections given in Table 2 are 
the sum of: (1) all the incremental corrections applied during the iterative procedure, (2) the 
corrections applied to initially align the series with each other, (3) the corrections applied to place 
each series within that particular IERS reference frame defined by the IERS combined Earth 
orientation series EOP(IERS) 08 C 04, and (4) the additional, common correction applied to the 
combined series in order to align the combined series with the IERS reference frame. The values 
of the uncertainty scale factors given in Table 2 are the products of all the incremental scale factors 
determined during the iterative procedure. The uncertainties of the bias-rate corrections given in 
Table 2 are the la standard errors of the incremental bias-rate corrections determined during the 
last iteration. There are no bias-rate entries in Table 2 for components that were either not used 
(e.g., the UTCSR SLR UT1 component) or not available (e.g., the NOAA IRIS Intensive polar 
motion components). Note that the same common IERS rate correction is applied to all the data 
sets, including those such as the Westford-Mojave single baseline VLBI series, for which no 
relative rate correction could be determined. Therefore, the rate corrections given in Table 2 for 
those series for which no relative rate corrections could be determined are simply the common 
IERS rate correction, but given in the natural reference frame of that series. In these cases, 
uncertainties for the rate corrections are not given. Also note that the entries for the bias-rate 
corrections in Table 2 for the IGS combined series EOP(IGS) 95 P 02 and EOP(IGS) 00 P 03 are 
the same. Since these entries were initially given within the same IGS reference frame, they were 
merged, and a common bias-rate correction was determined for them. 

Combined EOP Series: SPACE2013 

A Kalman filter was used to combine the series listed in Table 1 after the bias-rate 
corrections and uncertainty scale factors listed in Table 2 had been applied to them. The resulting 
combined Earth orientation series, SPACE2013, consists of values (Figure 1) and la standard 
errors (Figure 2) for polar motion, Universal Time, and their rates spanning September 28, 1976, 
to June 30, 2014, at daily intervals; and it is available in versions for which the epochs are given 
at either midnight or noon. Leap seconds have been restored to UT1, and the effects of the solid 
Earth and ocean tides have been added back to the UT1 values using the same models for these 
effects that were originally used to remove them from the raw series: namely, the Yoder et al. 
(1981) model for the solid Earth tides and the Kantha et al. (1998) model for the ocean tides. 
However, semidiurnal and diurnal ocean tidal terms have not been added to and are therefore not 
included in the SPACE2013 UT1 values. 

Figure 3 is a plot of the difference between the SPACE2013 polar motion, UT1, and LOD 
values and those of the IERS combined Earth orientation series EOP(IERS) 08 C 04. These two 
series are very consistent with each other, especially after January 1, 2000 when the root-mean- 
square (nns) of their difference is only 0.051 milliarcseconds (mas) for the 
x-component of polar motion, 0.045 mas for the v-component, 0.019 milliseconds (ms) for UT1, 
and 0.018 ms for LOD. Prior to 2000, the difference between these two series exhibits greater 
variability and even some systematic behavior, particularly in the x-component of polar motion. 
This systematic behavior is due to differences in the approaches used here and by the IERS to 
correct the bias and rate of the individual series before they are combined. 
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UT1-UTC (ms) Polar Motion X (mas) 


A Combined Earth Orientation Series: SPACE2013 






Figure 1 . Plots of the x-component of polar motion (la), y-component of polar motion (lb), UT1-UTC (lc), and excess length-of-day (Id) as given by the 
combined Earth orientation series SPACE2013. The discontinuous changes in the plot of UT1-UTC are caused by the presence of leap seconds. Note that the 
UT1-UTC values displayed in (lc) include the tidal variations, whereas the excess length-of-day values shown in (Id) do not. 
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A Combined Earth Orientation Series: SPACE2013 
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Figure 2. Plots of the lcj formal errors of the x-component of polar motion (2a), ^-component of polar motion (2b), UT1-UTC (2c), and excess length-of-day (2d) 
as given by the combined Earth orientation series SPACE2013. The inserts within panels (2a), (2b), and (2c) show that component’s post-1984 uncertainties on 
an expanded scale with the same units: milliarcseconds (mas) for polar motion, milliseconds (ms) for UT1-UTC. 
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Figure 3. Plots of the difference between the IERS combined Earth orientation series EOP(IERS) 08 C 04 and SPACE2013 formed by subtracting the SPACE2013 
values from those of the IERS series. The difference between the .^component of polar motion is shown in (3a), the difference between the v-componcnt is 
shown in (3b), the difference between UT1-UTC is shown in (3c), and the difference between the excess length-of-day is shown in (3d). 







COMB2013 


COMB2013 extends SPACE2013 by additionally incorporating the optical astrometric 
polar motion and UT1 series that was determined at the Bureau International de l’Heure (BIH) 
from an analysis of time and latitude observations by Li (1985; also see Li and Feissel 1986). This 
BIH optical astrometric series consists of values and uncertainties for polar motion and UT 1 that 
span from January 5.0, 1962 to December 31.0, 1981 at 5-day intervals. 

Data Preprocessing and Treatment of Tide-Induced Rotational Variations 

The BIH optical astrometric series was first preprocessed by removing leap seconds from 
the UT1 values and by correcting the UT1 values to be consistent with the extended definition of 
GST, as adopted by the IERS (IERS 1997, p. 149). Rotational variations caused by solid Earth 
and ocean tides were also removed from the UT1 values. The same models that were used to 
remove the tidal effects from the series combined to form SPACE2013 were also used to remove 
them from the BIH series: namely, the Yoder et al. (1981) model for the solid Earth tides and 
the Kantha et al. ( 1 998) model for the Mf, Mf, and Mm ocean tides. However, since the BIH UT 1 
measurements represent an average value over a 5 -day-long observation window, and since 5 days 
is a substantial fraction of the monthly and shorter-period tides, the amplitudes of these solid Earth 
and ocean tidal terms were attenuated prior to their removal from the BIH UT1 measurements. 
(See Gross 1996, p. 8735 and Gross et al. 1998, pp. 226-227 for further discussion about this 
point.) 

Adjustments Made to BIH Series Prior to Combination 

The preprocessed BIH optical astrometric series was combined with the space-geodetic 
series that comprise SPACE2013 after first: (1) correcting the BIH series to have the same bias, 
rate, annual terms, and semiannual terms as SPACE2013; (2) applying a constant multiplicative 
scale factor to the measurement uncertainties of the BIH series so that its residual, when 
differenced with SPACE2013, had a reduced chi-square of one; and (3) deleting those data points, 
if any, for which residual values were greater than four times their adjusted uncertainties. 
Unlike in previous years, software improvements have allowed the above adjustments 
to be detennined in a multivariate approach as was done for the series combined to form 
SPACE2013. The procedure used to determine these bias-rate and seasonal term corrections, 
uncertainty scale factors, and outlying data points has been described before (Gross 1996, pp. 
8735-8738) and will not be repeated here. The annual and semiannual terms of the BIH series 
were adjusted in order to correct for systematic, seasonally varying effects that are known to be 
present in optical astrometric measurements. Since the values of both the BIH series and the 
SPACE2013 series are given at midnight, interpolation error (see Gross et al. 1998, pp. 223-225) 
is automatically minimized when differencing these two series for the purpose of determining the 
adjustments to be made to the BIH series. Tables 3 and 4 give the resulting uncertainty scale factors 
and values and la standard errors of the corrections to the bias, rate, annual, and semiannual 
tenns thus detennined for the BIH series. When determining these uncertainty scale factors and 
the corrections to the bias, rate, and seasonal terms, no outlying data points were detected. 
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Combined EOP Series: COMB2013 


A Kalman filter was used to combine the BIH series with the adjusted space-geodetic series 
that comprise SPACE2013 after first applying to the BIH series the uncertainty scale factors and 
corrections to the bias, rate, annual, and semiannual terms that are given in Tables 3 and 4. The 
resulting combined Earth orientation series, COMB2013, consists of values (Figure 4) and la 
standard errors (Figure 5) for polar motion, Universal Time, and their rates that span from January 
20, 1962, to June 30, 2013, at daily intervals and is available in versions with epochs given at 
either midnight or noon. Leap seconds have been restored to UT1, and the effects of the solid 
Earth and ocean tides have been added back to the UT1 values using the same models for these 
effects that were originally used to remove them: namely, the Yoder et al. (1981) model for the 
solid Earth tides and the Kantha et al. (1998) model for the ocean tides. The full amplitude (i.e., 
no tidal terms attenuated) of the effects of the solid Earth and ocean tides at the epoch of the time 
tag were added back to the UT1 values. Semidiurnal and diurnal ocean tidal terms have not been 
added to and are, therefore, not included in the COMB2013 UT1 values. 


Table 3. Adjustments to Bias, Rate, and Stated Uncertainty of Optical Astrometric Series* 


Data Set 

Bias, mas 



Rate, mas/yr 


Uncertainty Scale Factor 

PMX 

PMY 

UT1 

PMX 

PMY 

UT1 

PMX 

PMY 

UT1 

BIH -2.456 ± 0.852 

-0.652 ± 0.699 

-7.541 ± 1.076 

1.050 ±0.526 

1.202 ±0.428 

5.583 ±0.673 

1.810 

1.677 

1.940 

ILS -50.509 ±2.356 

-3.212 ± 1.665 


0.427 ± 0.482 

-0.010 ±0.340 


2.158 

1.500 



*Reference date for bias-rate adjustment of BIH series is 1980.0. Reference date for bias-rate adjustment of ILS series is 1970.0. 


Table 4. Adjustment to Annual and Semiannual Terms of Optical Astrometric Series* 


Data Set 


Coefficient of Sine Term, mas 


Coefficient of Cosine Term, mas 


PMX PMY UT1 


PMX PMY UT1 


BIH 

Annual -5.426 ± 1.050 -7.390 ±0.882 5.185 ± 1.388 -2.869 ±1.121 9.711 ±0.922 -0.689 ±1.444 

Semiannual 2.292 ± 1.072 -0.537±0.897 -0.481 ± 1.415 0.872 ± 1.095 1.121 ±0.904 1.851 ± 1.415 


ILS 

8.072 ±3.323 -10.290 ± 2.344 

2.150 ± 3.312 1.811 ±2.335 


Annual -2.216 ± 3.287 8.236 ±2.315 

Semiannual 0.254 ± 3.296 8.771 ±2.322 


*Reference date for adjustment of annual and semiannual terms of BIH series is 1980.0. 
Reference date for adjustment of annual and semiannual terms of ILS series is 1970.0. 
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A Combined Earth Orientation Series: COMB2013 






Figure 4. Plots of the x-component of polar motion (4a), y-component of polar motion (4b), UT1-UTC (4c), and excess length-of-day (4d) as given by the 
combined Earth orientation series, COMB2013. The discontinuous changes in the plot of UT1-UTC are caused by the presence of leap seconds. Prior to the 
introduction of leap seconds in 1972, Coordinated Universal Time (UTC) was adjusted by introducing step and rate changes designed to keep it close to UT1 (e.g., 
IERS 1997, Table II-3), the effect of which is also readily apparent in (4c). Note that the UT1-UTC values displayed in (4c) include the tidal variations, whereas 
the excess length-of-day values shown in (4d) do not. 












POLE2013 


No optical astrometric observations taken at the stations of the International Latitude 
Service (ILS) were used when creating the BIH optical astrometric series used in COMB2013 (Li 
1985; Li and Feissel 1986). The ILS polar motion measurements (Yumi and Yokoyama 1980), 
which are based solely upon latitude observations made at the ILS stations are, therefore, 
independent of those comprising COMB2013 and have therefore been combined with them to 
form POLE2013. Being based solely upon latitude observations, the ILS series contains no UT1 
measurements, but consists solely of polar motion measurements that span 1899.8-1979.0 at 
monthly intervals. Although no uncertainties are given with the individual polar motion values, 
the precision with which they have been detennined is estimated to be 10-20 mas (Yumi and 
Yokoyama 1980, p. 27). An initial uncertainty of 15 mas was, therefore, assigned to each of the 
ILS polar motion values. Since this assigned measurement uncertainty will be adjusted later, its 
initial value is arbitrary as long as it is not zero and serves merely as an a priori estimate to be used 
in the series adjustment procedure described below. 

The ILS series was combined with COMB2013 to form POLE2013 after: 

(1) correcting the ILS series to have the same bias, rate, annual, and semiannual 
terms as COMB2013; 

(2) applying a constant multiplicative scale factor to the measurement uncertainties 
of the ILS series so that its residual, when differenced with COMB2013, had a 
reduced chi-square of one; and 

(3) deleting those data points, if any, whose residual values were greater than four 
times their adjusted uncertainties. 

Again, software improvements pennitted these adjustments to be detennined using a multivariate 
approach for fitting a bias, rate, and these seasonal terms to the difference of the ILS series with 
COMB2013 during 1962.0 to 1979.0. The measurement uncertainties of the ILS polar motion 
values were adjusted by detennining and applying a scale factor that made the residual of this fit 
have a reduced chi-square of one. During this procedure to determine uncertainty scale factors 
and bias, rate, and seasonal term conections, no outlying ILS data points were deleted since no 
data points had residual values greater than four times their adjusted uncertainties. Tables 3 and 4 
(in the COMB2013 section) also give the resulting uncertainty scale factors and values and la 
standard errors of the corrections to the bias, rate, annual, and semiannual terms thus determined 
for the ILS series. 

A Kalman filter was then used to combine the ILS series with the adjusted BIH and space- 
geodetic series that comprise COMB2013, after applying to the ILS series the uncertainty scale 
factors and corrections to the bias, rate, annual, and semiannual terms that are given in Tables 3 
and 4. The resulting combined Earth orientation series, POLE2013, consists of values (Figure 6a 
and 6b) and la standard errors (Figure 6c and 6d) for polar motion and its rate that span from 
January 20, 1900, to June 22, 2014, at 30.4375-day intervals. 
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A Combined Earth Orientation Series: POLE2013 
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Figure 6. Plots of the x-component of polar motion (6a), the y-component of polar motion (6b), the la formal errors of the x-component of polar motion (6c), and 
the la formal errors of the ^-component of polar motion (6d) as given by the combined polar motion series POLE2013. 






DISCUSSION 


The Kalman filter that was used here to combine Earth-orientation measurements contains 
a stochastic model of the process that is used to propagate the state vector and its associated state 
covariance matrix forward in time to the epoch of the next measurement. For polar motion 
excitation, the stochastic model includes a random-walk term with equal noise forcing both the x- 
and v-components of excitation and originally having a white-noise power spectral density of 
246.4 mas 2 /day (Morabito et al. 1988). This level of polar motion excitation process noise in the 
Kalman filter was increased to 739.2 mas 2 /day in order to better match the observed spectrum of 
polar motion excitation. Increasing the excitation process noise reduces the level of smoothing 
applied to the polar motion components of the propagated state vector and increases the covariance 
associated with those components. Thus, the SPACE2013, COMB2013, and POLE2013 polar 
motion and polar-motion excitation values are not as heavily smoothed as were those of 
SPACE2003, COMB2003, POLE2003, and earlier combinations produced at JPL; and the 
uncertainties assigned to the SPACE2013, COMB2013, and POLE2013 polar motion and polar- 
motion excitation values are somewhat larger. 

The LOD components of the IGS Final combined EOP(IGS) 95 P 02 and EOP (IGS) 00 P 
03 series were incorporated into JPL’s combinations starting with SPACE2007 and COMB2007. 
In general, GPS LOD estimates are contaminated by orbital artifacts that must be removed prior 
to combination with other LOD or UT1 measurements (Chin et al., 2009). When detennining the 
IGS Final combined LOD series, however, the IGS applies constraints that have the effect of 
removing these orbital artifacts. Therefore the IGS Final combined LOD series was incorporated 
into the Kalman filter as a true LOD measurement whose uncertainties are assumed to be “white” 
(i.e., random, uncorrelated in time, and with a flat power spectral density). Including the daily IGS 
Final combined LOD series in JPL’s combinations helps to compensate for less frequent UT1 
measurements and was found to improve the agreement of JPL’s combined LOD estimates with 
independent atmospheric and oceanic angular momentum values, especially at the highest 
frequencies. 

Since a Kalman filter has been used in generating SPACE2013, COMB2013, and 
POLE2013, the resulting polar motion and UT1 values are smoothed to a degree that depends upon 
both the spacing between the measurements being combined (which determines how far the state 
vector and state covariance matrix must be propagated forward in time) and the uncertainties that 
have been assigned to the measurements. Since improvements to the observing systems, both in 
the hardware and software and in the number of systems, have led to increasingly precise 
determinations of the Earth’s orientation, and since the time resolution of the measurements has 
generally increased in concert with the measurement precision, the degree of smoothing applied 
to the SPACE2013, COMB2013, and POLE2013 values is a function of time, with the earlier 
values being more heavily smoothed than the more recent values. 

Daily EOP values are reported in SPACE2013 since the NOAA IRIS and GSFC NEOS 
Intensive UT 1 values are given at daily intervals, as are the GPS and ILRS A combined SLR values 
(although gaps exist in each of these data sets). However, prior to the start of these data sets, the 
measurements combined to form SPACE2013 are given less frequently; therefore, the Kalman 
filter used to combine these measurements also interpolates them in order to produce a series of 
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equally spaced values. In order to be consistent with SPACE2013, daily EOP values are also 
reported in COMB2013 even though the BIH optical astrometric series used in COMB2013 is 
given at 5-day intervals. Thus, SPACE2013, COMB2013, and POLE2013 are equally spaced 
series of smoothed, interpolated Earth-orientation parameters. 

The combined Earth-orientation series SPACE2013, COMB2013, and POLE2013 are 
available from JPL’s Geodynamics and Space Geodesy Group via anonymous ftp: 
<ftp://keof.jpl.nasa.gov/keof/combinations/2013> 
and upon request from the authors: Todd.Ratcliff@jpl.nasa.gov or Richard.Gross@jpl.nasa.gov. 
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Acronyms and Terms 


ASI 

BIH 

CGS 

CODE 

COMB 


DSN 

EOP 

G3 

Gc 

GNSS 

GPS 

GSFC 

GST 

Hr 

HRAS 

ID 

IERS 

IGN 

IGS 

ILRS 

ILS 

IRIS 

ITRF 

IVS 

JPL 

Kk 

Ku 

LLR 

LOD 

mas 

MO 

ms 


Agenzia Spaziale Italiana 

Bureau International de THeure 

Centro di Geodesia Spaziale 

Center for Orbit Detennination in Europe 

extends the SPACE series by additionally incorporating BIH optical astrometric 
measurements of polar motion and UT 1 

Deep Space Network 

Earth orientation parameters 

IVS station code for 85-3 antenna at NRAO Green Bank, WV 

IVS station code for Gilmore Creek, Fairbanks, AK 

Global Navigation Satellite System 

Global Positioning System 

Goddard Space Flight Center 

Greenwich Sidereal Time 

IVS station code for HRAS 

Harvard Radio Astronomy Station at Ft. Davis, TX 
Identification 

International Earth Rotation and Reference Systems Service 

Institut National de lTnformation Geographique et Forestiere 

International Global Navigation Satellite Systems (GNSS) Service 

International Laser Ranging Service 

International Latitude Service 

International Radio Interferometric Surveying 

International Terrestrial Reference Frame 

International VLBI Service for Geodesy & Astrometry 

Jet Propulsion Laboratory 

IVS station code for Kokee Park, Kauai, HI 

IVS station code for 9-m antenna at Kokee Park, Kauai, HI 

lunar laser ranging 
length-of-day 

milliarcsecond 

IVS station code for 12-m antenna at Goldstone, CA 
millisecond 
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NASA 

NEOS 

NO A A 

NRAO 

NRCan 

Ny 

OCA 

PMX 

PMY 

POLE 


Ri 

rms 

SINEX 

SLR 

SPACE 

Sv 

T 

Ts 

USNO 

UT 

UTO 

UT1 

UTC 

UTCSR 

V 

VLB I 
VOL 

Wf 

Wz 


National Aeronautics and Space Administration 
National Earth Orientation Service 
National Oceanic and Atmospheric Administration 
National Radio Astronomy Observatory 
Natural Resources Canada 

IVS station code for Ny Alesund, Svalbard, Norway 

Observatoire de la Cote d’Azur 

polar motion X 
polar motion Y 

extends the COMB series by additionally incorporating ILS optical astrometric 
measurements of polar motion 

IVS station code for Richmond, FL 
root mean square 

Solution Independent Exchange 
satellite laser ranging 

a combination of independent space-geodetic measurements of the Earth's 
orientation 

IVS station code for Svetloe, Russia 

transverse component of Earth orientation from single-baseline VLBI 
IVS station code for 32-m antenna at Tsukuba, Japan 

United States Naval Observatory 
Universal Time 

Universal Time as determined at a particular observatory 

The principle form of Universal Time; mean solar time at 0° longitude obtained 

by correcting UTO for observatory location and polar motion. 

Coordinated Universal Time 

University of Texas Center for Space Research 

vertical component of Earth orientation from single-baseline VLBI 
very long baseline interferometry 
variation of latitude 

IVS station code for Westford, MA 
IVS station code for Wettzell, Germany 
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